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We propose a method to produce a definite number of ground-state atoms by adiabatic reduction 
of tlie depth of a potential well that confines a degenerate Bose gas with repulsive interactions. 
Using a variety of methods, we map out the maximum number of particles that can be supported 
by the well as a function of the well depth and interaction strength, covering the limiting case of 
a Tonks gas as well as the mean-field regime. We also estimate the time scales for adiabaticity 
and discuss the recent observation of atomic number squeezing (Chun et al., Phys. Rev. Lett. 95, 
260403 (2005)). 
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The controlled generation of many-body atomic num- 
ber states has been a long-standing goal in physics and 
success would open a door to a controlled study of en- 
tanglement 0, 0, 0| ; few-body tunneling and could 
also find important applications in quantum computing. 
One avenue towards this goal is the Mott insulator state 
where single-atom or multi-atom number states are pre- 
dicted [a, Q . Most experiments to date have used op- 
tical lattices where direct access to individual sites 
has not been accomplished and appears very difficult. A 
completely different approach was used in a recent ex- 
periment. The resulting atomic number squeezing was 
directly measured by atom counting, and a number state 
was inferred by accounting for known noise sources 0. 
In this Letter we analyze this new approach. 

The basic idea is to confine a degenerate Bose gas in 
an optical box with finite barrier height that can be con- 
trolled. The repulsive interaction between the atoms 
means that a finite box can only contain a maximum 
number of atoms. As the barrier height is slowly re- 
duced, atoms must leave, and the final number will be 
completely determined by the stopping point of the bar- 
rier. Since the confinement and the barrier in Ref. 
were realized by means of dipole optical traps we call 
this process "laser culling of atoms" although the same 
principle could also be implemented in other types of 
traps. The main theoretical questions that need to be 
addressed are how the maximum number depends on the 
potential parameter and interaction strength and how 
slow the potential should be changed in order to avoid 
excitations within the box. 

We address these questions by calculating and analyz- 
ing the energy levels of iV-particle ground states as well 
as excitations that are bound to the potential well. This 
is in general a very difficult task because of the interac- 
tions between the atoms. Therefore, we consider a simple 
model of one dimensional (ID) bosons with contact inter- 
actions in a square well potential. This model captures 
the essential features of the experimental system used 



recently to produce atomic number squeezing [9j. 

The letter is structured as following. First, we describe 
the model. Then we outline the idea of atomic culling in 
the limiting case of impenetrable bosons, which is fol- 
lowed by the discussion of another limit where modified 
mean- field picture is relevant. The diffusion Monte Carlo 
and direct diagonalization of the Hamiltonian bridge 
these two regions. Finally, we discuss the criteria for 
adiabaticity of the process and effect of the initial tem- 
perature. 

As usual, the interaction strength g, the coefficient 
of the (5-function interaction potential, depends on the 
three dimensional s-wave scattering length as well as 
the widths of the transverse wave functions a_L: g — 
2lT?as/ma\ when a_L S> Og [lo| . One can therefore vary 
g by changing the transverse confinement potential. The 
potential well depth Vb is controlled by the laser inten- 
sity for the barriers. The width L of the potential well 
is another parameter that can be adjusted over a wide 
region. Our model is made dimensionless by using the 
convention that the Planck constant h, atomic mass m, 
and the well width L are all unity. Thus, there are two 
dimensionless parameters in our model, the well depth in 
units of [h/Vf' jm and the interaction strength in units 
of h'^/mL. 

To explain how quantum culling happens we consider 
N impenetrable bosons in a ID well. In this situation the 
problem is mapped exactly to N non-interacting fermions 
in the same potential The eigenstates "^{xi, ■ ■ ■ Xn) 
of the problem are given by properly symmetrized set of 
single particle eigenstates (fij (x) 

*(a;i, • - -Xn) = S(f)j^{xi) ■ • •0j„(a;7v), (1) 

which in this situation for bosons is given by the abso- 
lute value of the Slater determinant. The ground state 
corresponds to indices j arranged from 1 to N. Since 
the finite ID well supports N bound single particle state 
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FIG. 1: Ground state energy of TV Tonks bosons in a square 
well. Zero of energy is at the top of the well. For Vb where 
two levels merge larger number of particles is not supported. 
These values are indicated by vertical dashed lines. 
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FIG. 2: Variational parameters /ti and hi2 that minimize the 
total energy for A'^ = 3 and g = 1. Dashed line shows the 
variational parameter for a wave function that uses only one 
variational parameter. 



only for 

Vo,N>{^{N-l)/Lf/2, (2) 

only then N impenetrable bosons can be trapped in it. 
As the well depth decreases to such a value that the most 
energetic single particle state delocalizes, only N ~1 par- 
ticles remain trapped in the well (see Fig.^. If initially 
only — 1 atoms were trapped, they remain trapped. 
Hence if initially there was uncertainty in the number of 
atoms in the trap [N oi N — 1) it is reduced. 

In the opposite limit of weak interaction, one may in- 
voke mean- field approximations. The simplest approach 
assumes that all the N particles are in the same single 
particle state, with its wave function determined vari- 
ationall y b y the standard Gross-Pitaevskii (GP) equa- 
tion |l2L [13 ■ This is only good for the ground state with 
all the particles deeply bound in the well. Near "ioniza- 
tion" threshold, only one particle should become weakly 
bound with long tails of the corresponding orbital reach- 
ing outside the well, while all the others remain tightly 
bound. This situation is better described by a Hartree- 
Fock wave function with — 1 particles in a state ipi 
and one particle in another state 4>2 0, ^3 ■ In order to 
minimize the total energy, (pi should be symmetric and 
nodelcss. For bosons, (t>2 docs not need to be orthogonal 
to 01, and should in fact also be symmetric and nodeless, 
because the A^-body wave function should be symmetric 
and nodeless with respect to any particle coordinate in 
the ground state. In Fig. |21 we show the result of vari- 
ational minimization of the total energy for the case of 
N = 3, and g = I, using variational functions of the form 

= v^exp(-Ki |a;|) (i = l,2). (3) 
We minimize the energy varying parameters ki and K2 



in the symmetrized wavefunction |^ for a given strength 
of interaction g and depth of the potential Vq . We allow 
only a single particle to populate the second orbital. Non- 
vanishing overlap between the orbitals must be taken into 
account. For a deep well Ki are comparable, which indi- 
cates a condensate state. As the depth of the well is 
reduced K2 tends to zero, which corresponds to delocal- 
ization of one particle. In comparison we show the result 
of minimization of total energy when only one variational 
parameter is used, then vanishing of the variational pa- 
rameter corresponds to delocalization of all the particles. 
For all the depths of the potential the energy is smaller 
when two-orbital wave function is used. 

In the limit of large number of particles and weak in- 
teraction the Thomas-Fermi approximation can be used. 
The approximation is valid when Ng ^ 1/L, where A^ is 
number of particles. Delocalization of the wave function 
in the standard GP equation in Thomas-Fermi limit can 
be expected when the chemical potential becomes equal 
to the level of the top of the well {fi = 0). This cor- 
responds to a depth of the potential Vb < g{N — 
When two asymmetric orbitals 4>i{x) and 1^2(3^) in the 
mean-field treatment |l4| are considered, minimization 
of energy functional leads to two coupled nonlinear equa- 
tions 

(hix) + giN - 2) + 2g l^aj") 0i = n^cfu , , 
/ 9\ (4) 

where h(x) is a single particle Hamiltonian. Neglecting 
kinetic energy term in the Hamiltonian and considering 
the situation when the orbital ^2 is almost delocalized, we 
see that the second equation is a single particle equation 
with potential modified by the rest of the atoms. As a 
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result N atoms are supported in the potential when 

Vn > 2g{N - 1)/L. (5) 

The energy calculated with two orbitals approach is 
smaller than calculated with a single orbital when Ng > 
LVq, which is consistent with the approximations. 

To describe the system for small values of g and small 
number of particles we resort to direct diagonalization 
of the many-body Hamiltonian. We limit the well by a 
box with infinite walls of size D ^ L. A single par- 
ticle wave functions i^j{x) can be either even or odd 
and inside of the well are proportional to cos{kx) and 
sin(A:a::) correspondingly. Outside of the well they are 
given by linear combination of exponents exp(±Kx), with 
K = \/k^ — 2Vo, and such coefficients that the wave func- 
tion is continuous at the boundary. When k is imag- 
inary the single particle states are extended. To find 
several single particle states, the values of k are adjusted 
so that the wave function vanishes at the boundary of 
the larger well. The N particle basis is constructed with 
these states. In the basis of second quantized states with 
Nj atoms on j's single particle state, single particle ener- 
gies Ej contribute to diagonal terms of the Hamiltonian 

H,=Y,Ejn„ (6) 

3 

while interaction gives contribution to other terms as well 
■^2 = ^ ^h^h^Ji^h UiM U \j3j4) , (7) 

jl ,32 ,33 ,34, 

{j.m\j3u)^ld.r,^rM^. (s) 

We diagonalize the full Hamiltonian H = Hi + H2- The 
bound states for N = 1,2,3 and g = 1 are shown in Fig.|2| 
for different Vq. 

As an alternative method to calculating the ground 
state of N bosons, we use the diffusion Monte Carlo ap- 
proach . Previously this method has been successfully 
applied to calculate the ground state of bosons with arbi- 
trary interaction. With this approach we also use a box 
of the size D ^ L. From the depths of the well Vq such 
that TV particles have smaller energy than iV — 1 particles 
we extrapolate smaller depth of the well and determine 
where the bound state of N particles is not bound any- 
more. Such a procedure leads to errors due to the finite 
time of the diffusion Monte Carlo evolution, finite size 
of larger box and extrapolation. In Fig. 0] we show re- 
gions where N particles can be supported. The errors are 
estimated to be comparable to the size of the symbols. 
Transition from iV = 1 to iV = 2, and from = 2 to 
A = 3 for 5 = 1 agrees quite well with diagonalization 
of the second quantized Hamiltonian. Transition from 




FIG. 3: Bound levels of bosons interacting with g — 1. Zero 
of energy is at the top of the well. Thick lines correspond 
to ground state of A'^ particles. Thin lines indicate bound 
excitation. Dashed-dotted line shows region when N — 2 
atoms are bound and A*' = 3 are not. Two vertical dashed 
lines indicate the regions where only A'^ = 1,2 particles are 
supported in the bound state (from left to right). Vertical 
dotted line shows the smallest Vb for which excited bound 
state is supported for = 3. 



A = 4 to A = 5 is close to what is predicted by Thomas- 
Fermi formula Eq. 

We now discuss how to achieve the atomic number 
states that we have studied above. We assume that our 
initial state is one with an unknown number of atoms but 
at essentially zero temperature. If the initial temperature 
is finite, a cooling procedure needs to be applied, which 
can be done by standard evaporation techniques. There- 
fore, one can prepare an initial state corresponding to 
one of the ground states that we have calculated above. 
However, the atomic number may not be the maximum 
allowed for the potential well, so that it is impossible 
to determine the initial atomic number from the depth 
of the potential well. One can then lower the barriers 
slowly until a threshold where it is no longer possible to 
hold the initial number of atoms. Then one atom will 
leave the well and there will be one less atom in the well. 
From this point on, as the potential barriers continue to 
reduce, there will be a one to one correspondence between 
the atomic number in the well and the intervals between 
the threshold values of the potential well. 

The discussion above assumes that initially system is 
not excited and that the excitations are not created dur- 
ing the process. There are two types of excitations: those 
when number of bound atoms is fixed and those when 
some atoms escape. For example, in Fig. |31 the tran- 
sition between points B and C corresponds to the first 
type when three bound atoms remain bound. On an- 
other hand, the transition indicated with arrows on the 
same figure corresponds to the second type when after 
excitation of two bound atoms one becomes unbound. 
All transitions between bold lines are of this type. When 
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one wants to control atoms with a single atom preci- 
sion these excitations are posing the main fundamental 
limitation. When bound state of A'' + 1 atoms becomes 
unsupported the state with N atoms has only one bound 
state. Hence, only unbounding excitation are possible at 
the final stage. To illustrate this with a specific example, 
examine Fig.OlandlSl There the level with = 2 atoms 
is highlighted with dashed-dotted line when = 3 has 
just become unsupported. Above it there are no excited 
states with N — 2. Also if one wants to obtain two atoms 
he better stop at the right side of the indicated inter- 
val, because on the right side the relevant excitation gap 
vanishes. First, we discuss the heating excitations (when 
number of atoms is preserved) in the connection with 
recent experiments, and later, consider the unbounding 
excitations during the final stage. 

During the evolution the potential barrier must vary 
sufficiently slowly in order to keep the system always in 
the ground state for each given number of atoms. The 
Planck constant divided by the excitation energy gaps 
give a time scale during which the wave function cannot 
change significantly for the evolution to remain adiabatic. 
For this problem the wave function undergoes significant 
changes in the range of the potential from one atom de- 
localization to another. 

In the recently reported experiment (21] j sub-Poissonian 
number statistics were directly observed by atom count- 
ing. The measured variances were nearly a factor of two 
below the shot-noise and could be attributed to known 
sources of technical noise. These observations were there- 
fore consistent with the production of number states. 
The change in atom number uncertainty was studied for 
different rates with which the depth of the potential is 
changed. Increase in uncertainty for larger rates can 
be attributed to generation of bounded excitation. We 
took the shape of the potential change as in the paper, 
Vdcpth(0 = Vbe"*/"^, and assumed the dependence of the 
number of atoms on the depth of the potential to be in the 
form A^ cx Vj^^pthJ with power a of order of one. We use 
GP equation and estimate the time scale of the excita- 
tion tgap from the energy of the Bogoliubov excitation of 
the longest wavelength [Tg] uj — (pUo/m)^^'^?:/ L^, where 
Uq — ATrha/m, a is the scattering length, m is the atom 
mass, p is the density, extension of the wave function 
in the non-confined direction. We would like to empha- 
size that the numbers below may be considered only as 
an order of magnitude estimate since as discussed above 
the delocalization depth may be different by some factor 
from what expected with GP equation. We estimated 
the average density by taking the parameters of the trap 
and numerically finding the ground state of the 3D GP 
equation, changing the degree of nonlinearity to assure 
confinement in the well. For the parameters of the ex- 
periment for A'finai we find tgap ~ 5 ms, and the evolution 
is adiabatic when t 3> 60 ms, consistent within an order 
of magnitude observed in the experiment. If one changes 



the well faster, excitations to the levels corresponding at 
the end to smaller number of atoms are unavoidable, re- 
sulting in larger spread in the distribution of number of 
atoms after the potential change as observed experimen- 
tally. 

As discussed above, the unbounding excitations put 
the most fundamental restriction when control with pre- 
cision of one atom is required. We consider unbounding 
excitations during the final stage in two limits. In the 
limit of strong interaction (Tonks gas), the system be- 
haves like a non-interacting fermion gas, and the excita- 
tions correspond to promoting the atoms from the filled 
single particle levels to the empty ones. In Fig. [3 we plot 
the energies of the lowest excited state together with the 
ground state energy for each A". During the final stage 
the relevant excitation energy gap is energy difference be- 
tween A^- and (A^— l)-particle states when A^-|- 1 particle 
have just stopped being supported. We illustrate this gap 
for N — 2 particles on Fig.jS] It can be calculated as the 
difference of two solutions of a transcendental equation. 
In the limit of large N it is £;gap ^ Tr'^(N+l/2). We would 
like to note that for A^ = 1 , 2 the similar linear depen- 
dence -Egap = AN + B holds except the slope A is approx- 
imately two times smaller. The depth of the well below 
which A^ particle cannot be bound is given by Eq. |21 As 
discussed above, the wave function varies significantly in 
the range of depth yo.w+i-Vb.AT = tt'^{N -1/2) / L'^ . The 
condition on potential depth rate of change, r = VQ{t)^ 
to be adiabatic is then (Planck constant here is unity as 
above) 

TT— < Vo^N+i - Vb,Ar. (9) 

-figap 

It follows that for large A^ the maximum rate is propor- 
tional to A^^ This can be understood from the Fig. |S1 
both gaps and gap intervals of Vb,7v are proportional to 
A^. As a result is in this limit when one goes to smaller 
and smaller number of atoms it becomes more and more 
difficult to remain adiabatic. In the case of A^ = I the 
maximum rate (in dimensionless units) is approximately 
unity (it does not vanish, because both gap and the depth 
intervals do not vanish). 

In the opposite limit of weak interaction, where the 
mean field picture is applicable, the energy gap is given 
by the energy difference between A^ and A^ -|- 1 particles 
when the latter stops being supported. For large A^ it 
is approximately -Egap = 3gAf^/2. The range of Vq for 
which A^ particles are supported in this limit is A^ in- 
dependent Vo,Ar+i — Vb,Af — 2g/L (see Eq. (SJ. Hence in 
this limit maximum rate r cx g^N"^ . It also becomes more 
and more challenging to obtain smaller number of atoms. 
The quadratic dependence on g shows that it is advan- 
tageous for the experimentalists to increase the effective 
interaction to make the maximum allowed rate larger. 

Finally, we would like to mention effect of initial exci- 
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FIG. 4: Parameter space for A'^ atoms to be bound by the 
well. Regions from bottom to top A'^ = 1, 2, 3, 4, and N > 5. 
Inset shows change in number of atoms along vertical dashed 
line. The squares are obtained with diffusion Monte Carlo 
method, the lines are guides for an eye. Two crosses for g = 1 
are obtained with diagonalization of full Hamiltonian. The 
inclined dashed line is the transition between A'' = 4 and 
N = 5 atoms in Thomas- Fermi limit, Eq. 0. 



tation for the process. As seen of Fig. O and Fig. [3 one 
may start from thin lines corresponding to excitations in 
TV > 1 systems and read the final state with = 1. In 
general, in both limits, if the final goal is to get N bound 
atoms and one starts with state with M atoms, M — N 
lowest many-body bound excitations are allowed. This 
means that for M ^ N one does not even have to start 
in degeneracy limit for the process to be possible. 

In conclusion, we have considered a system of N inter- 
acting bosons in a finite well. In such system a standard 
mean-field approach is not applicable. We have used dif- 
fusion Monte Carlo approach to obtain the parameters 
for N atoms to remain bound. The calculations agree 
with direct diagonalization for small number of atoms 
and analytical formulas in the limiting cases. As the 
depth of the well changes too fast excitation to levels 
that evolve into smaller number of particles can occur. 
We have estimated the critical rate for the recent exper- 
iments, which agreed within an order of magnitude. The 
limitations due to excitation during the last stage when 
precision within a single atom is required are discussed. 
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dation and the US Office of Naval Research, Quantum 
Optics Initiative, Grant N0014-04-1-0336. 
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FIG. 5: Bound levels of impenetrable bosons in a square well 
in dependence on the depth of the well. Zero of energy is at 
the top of the well. Thick lines correspond to ground state 
of N particles (same as Fig. . Thin lines indicate excita- 
tions. Dashed-dotted line shows region when N = 2 atoms 
are bound and A*' = 3 are not. Two arrows show the gap 
relevant during the final stage. If system starts at points A, 
B, or C it adiabatically evolves to point E {N = 2). Starting 
at point D it evolves to point F (A' = 1). A dashed level 
indicates excited state of A = 2, 3, and 4 configuration above 
which the system adiabatically evolves into A = 1 state. Ver- 
tical dashed lines show the regions where only N = 1,2,3 
particles are supported (from left to right). 
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